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Abstract 


Water is one of the critical needs of human life and living things. Therefore, proper planning is important 
for its consumption all over the world, especially in Iran. This study was conducted to evaluate the impacts 
of land use change on water yield in the Teraz watershed, Khuzestan province, Iran, from 1990 to 2050. 
Water yield was evaluated using integrated valuation of ecosystem services and tradeoffs (InVEST) tool 
and land use was mapped using the maximum likelihood classification in the ENVI 5.3, and CA-Markov 
in the TerrSet environment. Moreover, the relationship between water yield and landscape metrics, 
including the number of patches (NP), patch density (PD), landscape shape index (LSI), and Largest Patch 
Index (LPI), was considered based on the Geographically Weighted Regression (GWR) method in the Arc 
GIS 10.5. The results showed that the area of forest and rangeland in the studied area decreased in the last 
30 years, by 3199 and 1611 ha, respectively, and the area of agriculture and construction land uses increased 
by 4388 and 387 ha, respectively. It is predicted that in the next 30 years, 2442 ha of forests will decrease 
and the area of agriculture, rangeland, and construction land uses will increase by 1651, 687, and 102 ha, 
respectively. It was found that the total volume of available water yield in the region is equal to 26.5 Mm? 
in 2020, on average of 857 m° ha", and based on that, the ranking of land uses is as follows: 1. Construction, 
2. Rangeland, 3. Forest, 4. agriculture. The results of GWR confirmed that water yield had a significant and 
negative spatial relationships with the NP, PD, and LSI metrics (R?>=0.83, p-value>0.05), while its 
relationship with the LPI metric was a significant and positive relationship (R?>0.84, p-value>0.05). It can 
be concluded that the landscape features and land use pattern can determine the production and yield of 
water in the study area. 
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حکیده 


آب یکی از نیازهای اساسی حیات انسان‌ها و موجودات زنده به شمار می‌رود و برنامه‌ریزی مناسب برای 
استفاده از آن در سراسر جهان» بخصوص ایران. حائز اهمیت است. اين مطالعه با هدف ارزیابی اثرات 
کاربری اراضی بر خدمات اکوسیستمی تولید آب در حوضه آبخیز تراز» استان خوزستان در بازه زمانی 
۰ تا ۲۰۵۰ انجام شد. در اين راستاء ارزیابی تولید آب با استفاده از INVEST‏ و نقشه‌سازی کاربری 
اراضی با استفاده از روش طبقه‌بندی حداکثر احتمال در محیط 5.3 ENVI‏ و اتوماتای سلولى زنجیره 
مارکف در محيط 1617561 صورت پذیرفت. همچنین رابطه تولید آب با سنجه‌های سیمای سرزمین شامل 
تعداد. تراک شکل و مساحت لکه بر مبنای روش رگرسیون وزنی مکانی (GWR)‏ در محیط Arc GIS‏ 
5 مورد ارزیابی قرار گرفت. نتايج نشان داد در ۳۰ سال گذشته مساحت جنگل‌ها و مراتع در منطقه 
مطالعاتی به ترتیب حدود ۳۱۹۹ و ۱2۱۱ هکتار کاهش و مساحت اراضی کشاورزی و ساخته شده‌ها به 
ترتیب ۳۸۸ و ۳۸۷ هکتار افزايش یافته‌لند. پیش‌بینی می‌شود که در ۲۰ سال آینده نیز ۲۶۶۲ هکتار از 


مساحت جنگل‌های منطقه کاسته شود و مساحت کاربری‌های کشاورزی, مراتع و ساخته شده‌ها به ترتیب 
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۱ ۸۷ و ۱۰۲ هكتار افزايش يابد. كل حجم آب توليدى قابلدسترس در منطقه در سال ۲۰۲۰ معادل 
۵ میلیون مترمکعب (به‌طور متوسط AOV‏ مترمکعب در هکتار) و میزان تولید آب در ساخته‌شده‌هاء مراتع» 
جنگل‌ها و اراضی کشاورزی به ترتیب ۰۳۱۲۳ ۰۱۰۷٤‏ ۱۰۷۱ و ۳۷۶ مترمکعب در هکتار برآورد گردید. نتایج 
GWR‏ تأييد کرد که تولید آب با سنجه‌های تعداد لکه» تراکم لکه و شاحص شکل سیمای سرزمین رابطه 
مکانی معنادار و معکوس دارد (R20.83, p-value>0.05)‏ درحالی‌که رابطه Lol‏ شاخص بزرگ‌ترین لکه 
يك رابطه معنادار و مستقیم است p-value>0.05)‏ ,87-0.84). بر اين اساس می‌توان كفت که ویژگی‌های 
ساختاری سیمای سرزمین و الگوی کاربری اراضی حاکم بر منطقه می‌تواند تعیین‌کننده تولید آب باشد. 
کلیدواژه‌ها: خدمات اکوسیستمی. تولید آب. تغییرات کاربری اراضی» اینوست. رگرسیون وزنی مکانی. 

۱- مقدمه 

وضعیت و روند تغییرات خدمات اکوسیستمی در جهان. حیات انسان و موجودات دیگر را به‌طور مستقیم و 
غيرمستقيم تحت تأثير قرار می‌دهد و تهدیدی جدی برای امنیت محیطزیستی آینده و رفاه جوامع کنونی به شسمار 
می‌رود. خدمات اکوسیستمی تحت تأثير بحران‌های محیط‌زیستی شامل تغییرات اقلیمی و کاربری اراضی. رشد 
جمعیت. توسعه سریع اقتصادی و بهره‌برداری بی‌رویه از منابع به‌شدت در حال کاهش هستند؛ آزاین‌رو نیاز است که 
شرایط و وضعیت خدمات اکوسیستمی به‌طور مناسب و دقیق درک شود تا بر اساس أن بتوان رویکردهای مناسب 
مدیریتی در راستای استفاده مناسب و حفاظت خدمات اکوسیستمی اتخاذ کرد (احمدی میرقائد ۱۳۹۷). 

فعالیت‌های اقتصادی - اجتماعی انسان تغییرات اساسی در ساختار و عملکرد اکوسیستم‌ها ایجاد می‌کند که در 
نتيجه آن خدمات اکوسیستمی مختلف تحت تأثير قرار می‌گیرند. توسعه کاربری‌ها و تغييرات آن‌ها ازجمله مهم‌ترین 
اقدامات انسانی است که ویژگی سیستم‌های اکولوژیک را تغییر می‌دهند و تأثیرات مهمی بر خدمات اکوسیستمی 
ایجاد می‌کنند (احمدی ميرقائد و SOL See‏ ۲۰۲۰). تغییرات کاربری اراضی به‌عنوان نماینده تغییرات حدمات 
اکوسیستم برای ارزیابی توسعه و تغییرات محیط‌زیستی آینده در نظر گرفته می‌شسود. ازاین‌رو توجه به تغییرات 
کاربری اراضی و مدیریت آن‌ها حائز اهمیت است SY)‏ و سانگ"؛ ۲۰۱۹؛ پنگ و همکاران" ۲۰۱۹). 

لی و همکاران* (۲۰۲۱) اثرات تغییرات اقلیمی و کاربری بر خدمات اکوسیستمی در حوضه ایلی بلخاش در 
آسیای مرکزی را موردمطالعه قرار دادند. آن‌ها تغییرات مکانی-زمانی در آب‌وهواء کاربری زمين و چهار حدمات 
اکوسیستمی شامل تولید آب» حفاظت خاک ذخیره كربن و کیفیت زیستگاه‌ها را از ۱۹۹۲ تا ۲۰۱۶ ارزیابی کردند. 


مونپراپوسورن* (۲۰۱۷) چگونگی پیش‌بینی خدمات آینده اکوسیستم ساحلی در استان ساموتساکورن را بر اساس 
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تغییرات كاربرى و اقليمى موردمطالعه قرار داد. وى نشان داد كه افزايش سطح درياء دما و بارش با کاهش پوشش 
جنگلی و درنهايت تضعيف عملكرد و خدمات اكوسيستم ساحلى همراه است و از تاب‌آوری اكوسيستمها جلوكيرى 
می‌کند. شيريكى و همكاران' (۲۰۱۷) تأثيرات آتى تغييرات اقليمى و كاربرى اراضى را بر خدمات اكوسيستم 
جندمنظوره در دره استوباى اتريش در آلب مركزى را مورد ارزيابى قرار دادند. آن‌ها شش خدمات اكوسيستمى شامل 
توليد علوفه. كيفيت علوفه. حاص‌لخیزی خاک كيفيت آب. ذخيره كربن و ارزش‌های زيباشناختى را بر اساس 
مدل‌های مبتنى بر ویژگی‌های گیاهی برای شرايط فعلى و آينده (در سالهاى ۲۰۵۰ و ۲۱۰۰) وبا در نظر گرفتن 
شرايط اقتصادى-اجتماعى موردبررسی قرار دادند. هی و همکاران" (۲۰۱۷) اثرات تغيير کاربری بر كيفيت زيستكاه را 
با استفاده از چارچوب یکپارچه شبیه‌سازی اتوماتای سلولی و ابزار اینوست" در منطقه وانژو واقع در جين انجام دادند. 
آن‌ها کیفیت زیستگاه‌ها را طی سه سناریوی گسترش سریع شهری. توسعه شهری هوشمند و حفاظت محیطزیست 
بررسی کردند. وو و همکاران" (۲۰۲۰) خدمات اکوسیستمی وابسته به آب شامل عملکرد آب. تولید رسوب و تصفیه 
آب را در تعامل با تغییرات کاربری اراضى و اقلیمی در حوضه رودخانه لوانهه جين با استفاده از ابزار اينوست 
موردمطالعه قرار دادند. پلنگ و همکاران" (۲۰۲۱) خدمات تأمين آب در حوضه رودخلنه زرد واقع در جين را بر 
اساس مدل اینوست و روش تحلیل سناریو انجام دادند. آن‌ها از داده‌های پوشش و کاربری زمین» هواشناسی و نوع 
Se‏ به‌عنوان ورودی برای تجزیه‌وتحلیل الگوی مکانی و زمانی عملکرد آب و تأثیر بارش و تغییر کاربری زمين بر 
تولید آب استفاده کردند. اقصائی و همکاران" (۲۰۲۰) OLA‏ تغییرات کاربری و پوشش زمين را بر هیدرولوژی 
حوضه و عملکرد رسوب در حوضه آبریز انزلی» ايران را موردمطالعه قرار دادند. آن‌ها از مدل ارزیابی خاک و آب 
سوآت" برای ارزیابی تأثیرات تبخیر و تعرق» تولید آب و تولید رسوب استفاده کردند. 

حوضه آبخیز رودخانه تراز واقع در استان حوزستان در بخش نیمه‌عشک ايران واقع شده است و اثرات انسانی را 
می‌توان در تغيبر ویژگی‌های طبیعی آن بدخوبى مشاهده کرد. در بخش‌های شمالی منطقه میزان توسعه ساخته‌شده‌ها 
نسبت به بخش‌های جنوبی کمتر اتفاق افتاده است. در بخش‌های جنوبی و غرب منطقه choy‏ انسان و توسعه بسیار 
چشمگیرتر است و انجام فعالیت‌های انسانی ازجمله کشاورزی و توسعه زیرساخت‌ها در اين مناطق بسیار شدید 
است که باعث تخریب محیط‌زیست و کاهش کیفیت خدمات اکوسیستمی آن شده است. ازاينروء مطالعه اثرات 


1 Schirpke 

2 He et al. 

3 Integrated Valuation of Ecosystem Services and Tradeoffs (INVEST) 
4 Wu et al. 

5 Yang et al. 
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7 Soil and Water Assessment Tool (SWAT) 
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منابع طبیعی آن حائز اهمیت است. اين مطالعه با هدف (۱) مدل‌سازی خدمات اکوسیستمی تولید آب در حوضه 
مطالعاتی دو بازه زمانی ۱۹۹۰ 6 ۲۰۵۰ و (۳) ارزیایی ارتباط بين تولید آب با تغییرانت سیمای سرزعین با اسفاده از 


سنجه‌های سیمای سرزمین صورت پذیرفت. 


۲- مواد و روش‌ها 

۱-۲- منطقه مطالعاتى 

حوضه آبخیز رودخانه تراز با مساحتی بالغبر ۳۲۵۰۰ هکتار در شمال شرقی استان خوزستان در ايران قرار دارد 
(شکل ۱). اين حوضه یکی از زیر حوضه‌های اصلی رودخانه شور است که درنهایت به سد گتوند منتهی می‌شود. 
میانگین بارندگی و دمای سالیانه در حوضه به ترتیب WE‏ میلی‌متر و ۲۳ درجه سانتی‌گراد برآورد شده است. حوضه 
آبخیز تراز از نظر توپوگرافی متنوع است به‌طوری‌که در قسمت‌های شمالی پستی و بلندی‌های شدید و در نواحی 
جنوبی و غربی مناطق نسبتاً هموار دیده می‌شود. جنگل‌ها بخش‌های شمالی و شرقی حوضه آبخیز را در برگرفته‌اند 
درحالی که کشاورزی و مرتع کاربری‌های عمده در بخش‌های جنوبی به شمار می‌روند. توپوگرافی و پوشش ALS‏ 
متنوع باعث افزایش تنوع اکولوژیکی در منطقه شده است كه در نتيجه آن فراهم‌سازی خدمات اکوسیستمی مختلفی 
را به دنبال دارد. بااينحال» توسعه نامناسسب. بهره‌برداری زياد از منابع و مدیریت ضعیف daol‏ عمدتاً در قسمت 


جنوبی منطقه موردمطالعه خدمات اکوسیستمی به‌ویژه خدمات مرتبط با آب را به‌شدت تحت تأثیر قرار داده است. 
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شکل ۱- موقعيت جغرافیایی منطقه مطالعاتی 
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-Y-Y‏ جمع‌آوری داده‌ها 

ليست داده‌های مورداستفاده در اين مطالعه در جدول ۱ نشان داده شده است. نقشه‌های کاربری اراضی بر مبنای 
داده‌های چند طیفی ماهواره لندست ۸ و ۵ و با استفاده از روش نظارت شده حداکثر احتمال" تولید شدند. آمار 
ایستگاه‌های کلیماتولوژی منطقه ایجاد شده توسط سازمان هواشناسی ايران برای تهیه داده‌های اقلیمی موردتوجه قرار 
گرفت. زیرحوضهها نيز بر اساس مدل رقومی ارتفاعی" وبا استفاده از ابزار ارك هيدرو" در 10.5 Arc GIS‏ 


نقشه‌سازی شدند. 


حدول ۱- داده‌های مورداستفاده در اين مطالعه و توضیحات آن‌ها 


معیار واحد نحوه تهیه و آماده‌سازی 
مدل رقومى ارتفاع m (DEM)‏ استخراج شده از سایت سازمان زمین‌شناسی ایالاتمتحده امریکا CUSGS*)‏ 
توليد شده بر مبنای تصاوير لندست ۸ و ۵ و با استفاده از روش حداكثر احتمال در 
کاربری اراضی 
محيط نرم‌افزار5.3 ENVI‏ 
ميانكين بارندگی ساليانه mm‏ استخراج شده از آمار ایستگاه‌های هواشناسی منطقه (سازمان هواشناسی کشور) 
ميانكين حداكثر و حداقل دمای روزانه °C‏ استخراج شده از آمار ایستگاه‌های هواشناسى منطقه (سازمان هواشناسی كشور) 
درصد رس» سيلت و ماسه خاک i‏ روش هيدرومترى 
رع OC %OM = %OC x 1.72 i EE‏ كربن آلى خاک است (درخشان بابايى و همکاران 
۳۲۰( 
تبخير و تعرق مرجع mm (ETo)‏ محاسبه شده بر اساس معادله ٤‏ (هارگریوز OAIE‏ 
عمق لايه محدودکننده mm‏ نقشه عمق خاک منطقه مطالعاتی 
ضريب آب قابلدسترس كياهى mm (PAWF’)‏ محاسبه شده بر اساس معادله ۷ SY)‏ و همکاران“ ۲۰۱۷) 
ضريب تبخير و تعرق كياهى (Ke)‏ استخراج از منابع علمى (شارپ و همکاران؟» ۲۰۲۰) 
عمق ريشهيابى استخراج از abe‏ علمى (شارپ و همکاران ۲۰۲۰) 
فاکتور فصلى mm (Z)‏ محاسبه شده بر اساس وقايع بارندكى منطقه در طول يك سال 
نقشه زیرحوضه‌ها استخراج شده بر اس DEM‏ 52 محيط Arc Hydro‏ 


1 Maximum Likelihood Classification 
2 Digital Elevation Model (DEM) 

3 Arc Hydro 

4 United States Geological Survey 

5 Derakhshan-Babaei et al. 

6 Hargreaves 

7 Plant Available Water Fraction 

8 Lang et al. 

9 Sharp et al. 


۳ جغرافيا و مخاطرات محيطى شمارة دوم 


۳-۲- ارزیابی خدمات اکوسیستمی تولید آب! 
در مورد چگونگی تغییرات الگوهای کاربری اراضی بر تولید سالانه آب‌های سطحی ارائه می‌دهد. اين مدل بر اساس 
می‌کند. بعلاوه» میانگین و کل تولید آب را نیز در سطح زیرحوضه‌ها محاسبه می‌کند (شارپ و همکاران, ۲۰۲۰؛ 


احمدی میرقائد. ۱۳۹۷). 


AET x 
Px 


k= (1 - =) xk (\) 

كه در آن ,۸۳ و ,۴ به ترتيب ميزان تبخير و تعرق واقعى و بارندكى ساليانه را بر حسب ميلىمتر در پیکسل × 

نشان می‌دهند. برای هر نوع كاربرى اراضی» نسبت تبخير و تعرق در بيلان آبى بر اساس منحنى بوديكو است كه 
توسط فو و همکاران" (۱۹۸۱) و ژانگ و همکاران" (۲۰۰۶) ارائه شده است: 


AETy _ PETx _ PETx\"44/ 

3 85, OG) 0 
PET, = Ke(lx) x ETo (۳ 
ETo = 0.0013 x 0.408 x RA x (Tayg + 17) x (TD — 0.0123p)°76 (£) 
Wy = 7 2 + 1.25 5 


كه PET,‏ بیانگر تبخير و تعرق يتانسيل است و wx‏ یک پارامتر غير فيزيكى را نشان می‌دهد که مشخص کننده 
خصوصيات اقليمى و خاکی طبيعى است. ET‏ تبخير و تعرق مرجع است كه شرايط اقليم محلى را بر تبخير و تعرق 
پوشش كياهى مرجع مثل علف و يونجه را كه در آن محل رشد يافته بیان می‌کند. 166 ضريب تبخير و تعرق كياهى 
وابسته به يك نوع كاربرى (1) است كه تا حد زيادى توسط ویژگی‌های رويشى سرزمين و كاربرى اراضى تعيين 
ge‏ شود (آلن و همكاران*. ۱۹۹۸). RA‏ ميزان تابش فرازمينى بر حسب (MI Mm? d’)‏ است. Tavg‏ ميانكين حداكثر و 
حداقل دمای روزانه و TD‏ اختلاف ميانكين حداقل و حداكثر دماى روزانه بر حسب درجه سانتى كراد هستند. P‏ 
ميزان بارندگی ساليانه بر حسب میلی‌متر است. 2 يك ثابت تجربى است كه به‌عنوان "فاكتور فصلی" " در نظر كرفته 
مى شود و الكوى بارش محلى و ویژگی‌های هیدروژئولوژیکی را نشان می‌دهد. اين ثابت با تعداد رويدادهاى باران در 
سال N)‏ همبستگی مثبت دارد. حداقل مقدار 7 زمانى است که عمق ريشه صفر SL)‏ لخت) باشد و ply‏ با ۱/۲۵ 


1 Water yield 

2 Budyko 

3 Fuetal. 

4 Zhang et al. 

5 Allen et al. 

6 Seasonal factor 


سال دوازدهم اثر تغييرات كاربرى اراضى بر توليد اب در .... Yo‏ 


خواهد بود و مقدار حداكثر آن به ۵ محدود مى شود (دناهو و همکاران؛ ۲۰۱۲؛ یانگ و AWC? .)۲۰۰۸ Tab Rae‏ 
محتوای حجمی آب قابل‌دسترس گیاهی است که بر اساس بافت خاک و عمق ریشه‌زایی مؤثر تعيبن می‌شود. AWC‏ 
مقدار آب مورداستفاده توسط گیاه موجود در خاک را تعیین می‌کند و بر مبنای ظرفیت آب قابلدسترس گیاهی 
(PAWC')‏ و حداقل عمق لايه محدودکننده ريشه (Rid)‏ و عمق ریشه‌زایی (Rd)‏ پوشش گیاهی برآورد می‌شود. 
AWC, = Min(Rld, Rd) x PAWC a)‏ 
عمق لايه محدودکننده ريشه. عمق خاكى است که در آن نفوذ ريشه به دليل ویژگی‌های فيزيكى يا شيميايى مهار 
می‌شود. عمق ريشهزايى پوشش كياهى اغلب به‌عنوان عمقى در نظر كرفته مى شود كه در آن ۹۵ درصد از زیست‌توده 
ريشه كياهى وجود دارد. PAWC‏ ظرفيت آب قلبلدسترس كياهى است كه تفاوت بين ظرفيت مزرعه و نقطه 
يزمردكى را مشخص مىكند و ضريب آن PAWE)‏ بر اساس فرمول زیر برآورد مىشود: 


PAWF = 54.509 — 0.132 x SA% — 0.003 x (SA%)? — 0.055 x SI% — 0.006 x (SI%)? — 
0.738 x CL% + 0.007 x (CLM) — 2.688 x OM% + 0.501 x (OM%)? (VY) 


كه در SA ASI CLT‏ و ZOM‏ به ترتيب درصد محتواى رس سیلت» ماسه و مواد آلى را نشان می‌دهند 
Sb)‏ و همكاران', SY ١‏ و همکاران» ۷( 


4-7-اثرات تغيير کاربری اراضی بر تولید آب 
در اين مطالعه. يس از تهیه داده‌های موردنیاز (جدول ۱ با استفاده از مدل تولید آب در اینوست. اثرات تغییر 
کاربری اراضی در منطقه مطالعاتی در دو بازه زمانی ۲۰ ساله (۱۹۹۰ تا ۲۰۲۰ و ۲۰۲۰ تا ۲۰۵۰) بر تولید آب مورد 


ارزیابی قرار گرفت. برای اين منظورء نقشه‌های کاربری اراضی برای سال‌های ۱۹۹۰ و ۲۰۲۰ بر اساس تصاویر چند 
طیفی لندست ۵ و ۸ و با استفاده از روش حداکثر احتمال در محيط 5.3 ENVI‏ نقشه‌سازی شدند. همچنین با استفاده 
از روش اتوماتای سلولی - زنجیره مارکف" در محبط 1617566 نقشه کاربری اراضی برای سال ۲۰۵۰ مدل‌سازی 
گردید. طی سه سناریو تولید آب بر اساس نقشه‌های کاربری اراضی ۰۱۹۹۰ ۲۰۲۰ و ۲۰۵۰ محاسبه و میزان تغییرات 
تولید آب در دو بازه زمانی ۳۰ ساله برآورد گردید. لازم به ذکر است که دراین‌باره ساير ورودی‌ها در اجرای مدل 


ثابت در نظر گرفته شدند. 


1 Donohue et al. 

2 Yang et al. 

3 Available Water Content 

4 Plant Available Water Content 

5 Plant Available Water Fraction 

6 Yang et al. 

7 Lang et al. 

8 Cellular Automate — Markov (CA-Markov) 


۳۹1 جغرافيا و مخاطرات محيطى شمارة دوم 


0-۲- اثرات خردشدگی سیمای سرزمین بر تولید آب 

سنجه‌های سیمای سرزمین شاخص‌هایی هستند که جنبه‌های مختلف سیمای سرزمین را از نظر مساحت. شکل» 
ترکیب» توزیع و هندسه به‌صورت ارزش‌های کمی نشان می‌دهند که در نتیجه آن تحلیل» ارزیابی» يايش و طرح‌ریزی 
مکانی الگوی کاربری اراضی تسهیل می‌شود. بعلاوه» با کمک آن‌ها می‌توان میزان خردشدكى زیستگاه‌ها و سیمای 
سرزمین را ارزیابی کرد (احمدی میرقائد و همکاران؛ ۲۰۱۸؛ ينك و همکاران؛ VV‏ مک گاریگل و همکاران " 
۲ در این مطالعه. ارزیابی اثرات خردشدگی سیمای سرزمین بر تولید آب با استفاده از مهم‌ترین سنجه‌های 
سیمای سرزمین شامل تعداد لکه (NPE)‏ تراکم لکه (PD*)‏ بزرگ‌ترین لکه LPI)‏ و شکل سیمای سرزمین (LST)‏ 
مدنظر قرار كرفت که در جدول ۲ توصیف و تشریح شده‌اند. برای اين منظور» سنجه‌های سیمای سرزمین بر مبنای 
نقشه کاربری اراضی سال ۲۰۲۰ منطقه مطالعاتی و در محيط Fragstats‏ در سطح زیرحوضه‌ها محاسبه شدند و رابطه 


جدول -Y‏ سنجه‌های سیمای سرزمین و توصیفات آن‌ها Ke)‏ گاریگل و همکاران. ۲۰۰۲) 


سنجه واحد توصیف 
Ba‏ تعداد لکه‌های موجود در یک سیمای سرزمین مشخص را نشان می‌دهد. هر جه بزرگ‌تر از یک باشد 
00 خردشدكى بيشتر سيماى سرزمين را بیان می‌کند. 
تراكم لكه تعداد به ازای | تعداد لکه‌های موجود در سيماى سرزمين به نسبت كل مساحت آن را نشان می‌دهد. هر جه بزرگ‌تر از 
٠‏ هكتار يك باشد خردشدگی بيشتر سيماى سرزمين را بیان می‌کند. 
9 شاخصى است برای اندازهكيرى كل سيماى سرزمين بر حسب بزرگ‌ترین لكه و بين ١‏ تا ۱۰۰ متغير 
شاخص بزركترين لكه 1 ۳ p‏ 
است. هر جه بزركتر از يك باشد خردشدگی كمتر سيماى سرزمين را نشان می‌دهد. 
شاخص شكل سيماى شاخصى استاندارد برای اندازهكيرى كل لبه يا تراكم لبه در يك سيماى سرزمين مشخص است و هر 
سرزمين 0 جه بزرگ‌تر از يك باشد نشان‌دهنده پیچیدگی بيشتر مرز سيماى سرزمين است. 


1 Ahmadi Mirghaed et al. 

2 Peng et al. 

3 McGarigal et al. 

4 Number of Patch 

5 Patch Density 

6 Largest Patch Index 

7 Landscape Shape Index 

8 Geographically Weighted Regression 


سال دوازدهم اثر تغییرات کاربری اراضی بر تولید آب SO‏ بی ۳۷ 


۳- نتایج و بحث 

نقشه‌های کاربری اراضی منطقه برای سال‌های ۱۹٩۰‏ ۲۰۳۰ و ۲۰۵۰ در شکل ۲ مشخص شده‌اند. نتایج نشان داد 
که کاربری‌های جنگل و مرتع بیشترین مساحت نیمه شمالی حوضه را در برگرفته‌اند؛ درحالی که در نیمه جنوبی 
حوضه کاربری‌های کشاورزی, مرتع و ساخته شده‌ها گسترش یافته‌اند. کاربری‌های جنگل» مرتع» کشاورزی, ساخته 
شده‌ها و رودخانه‌ها به ترتیب ۱۰2۵۵ ۰۱۳۱۳۳ ۷۰6 ۶۳۲ و ۸۸۲ هکتار از مس احت منطقه را در سال ۲۰۲۰ در 
برگرفته‌اند. در سال ۱۹۹۰ مساحت اختصاصصيافته به کاربری‌های مذکور به ترتیب ۱۳۳۵۶ NEVER‏ ۰۳۲۱۶ 1۵ و 
۱ هکتار بوده است و بر اساس ادامه شرایط کنونی حاکم بر منطقه در ۳۰ سال آینده پیش‌بینی می‌شود که مقادیر 
مساحت‌های مذکور در سال ۲۰۵۰ به ترتیب به ۸۰۱۳ ۰۱۳۸۲۰ ۲۵۵ ۵۳۶ و ۸۸۷ هکتار تغییر یابند. بدين ترتیب 
مشخص شد که در ۲۰ سال گذشته مساحت جنگل‌ها و مراتع در منطقه مطالعاتی به ترتیب حدود ۳۱۹۹ و ۱۹۱۱ 
هکتار کاهش و مساحت اراضی کشاورزی و ساخته‌شده‌ها به ترتیب 1۳۸۸ و ۳۸۷ هکتار افزایش یافته‌اند. پیش‌بینی 
می‌شود که در ۳۰ سال آینده نيز ۲۶۶۲ هکتار از مساحت جنگل‌های منطقه کاسته شود و مساحت کاربری‌های 


کشاورزی. مراتع و ساخته‌شده‌ها به ترتيب AWN‏ ۷ و ۱۰۲ هکتار افزایش یابد. 


۳ 3 ۳ ع‎ ۳۳ te kı 
E كة ۰ رودخانه 88 انسان ساخت وس كشاورزى ا8ا مر تع 7۳۷ جتكل‎ 5 wo 


شكل ۲- نقشه‌های كاربرى اراضى منطقه مطالعاتى در سالهاى 1۹۹۰ ۲۰۲۰ و ۲۰۵۰ 


شکل ۳ روند تغییر مساحت اختصاصيافته به کاربری‌ها در منطقه مطالعاتی در سال‌های ۰۱۹۹۰ ۲۰۲۰ و ۲۰۵۰ 
را نشان می‌دهد. نتایج نشان می‌دهند که طی ۳۰ ساله گذشته توسعه و فعالیت‌های کشاورزی در منطقه با شدت زياد 
و بدون توجه به ظرفیت اکولوژیک و شرایط محیط‌زیستی منطقه انجام شده است که در بى آن اراضی جنگلی و 
مراتع تخریب و به کشاورزی و ساخته‌شده‌ها تبدیل شده‌اند. اين عوامل باعث افزایش خردشدگی زیستگاه‌ها و از 
هم گسیختگی سیماهای سرزمین حوضه شده است. درصورتی که همین روند و شرایط در توسعه کاربری‌ها در منطقه 
ادامه داشته باشد پیش‌بینی می‌شود که در ۳۰ سال آینده نيز بيش از ۲۵ درصد مساحت اراضی جنگلی تخریب و به 


کاربری‌های bas So‏ شود. شکل ؛ نيز روند تغییرات کاربری اراضی بر مبنای سنجه‌های سیمای سرزمین شامل 


YA‏ جغرافيا و مخاطرات محيطى شمارة دوم 


تعداد لکه» تراكم لکه» شاخص شكل سيماى سرزمين و شاخص بزركترين لكه را در بازه زمانى ۳۰ ساله (۱۹۹۰ تا 
۰ در مقياس زیرحوضه‌ها مشخص کرده است. نتايج حاكى از آن است که در ۳۰ سال اخير تعداد لکه تراكم 
لكه و شاخص شكل سيماى سرزمين در زیرحوضه‌ها روند افزايشى داشته درحالىكه شاخص بزرگ‌ترین لكه کاهش 
يافته است. اين موارد نشان می‌دهد كه خردشدكى و كسستكى لکه‌های كاربرى در منطقه طى ۳۰ سال اخیر افزايش 
يافته است که ناشى از توسعه فعاليتهاى كشاورزى و افزايش ساخته‌سازها بخصوص توسعه جاده و راه است. 
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ساخته شده ها کشاورزی 


شکل ۳- مساحت اختصاص يافته به هر کاربری در منطقه مطالعاتی در سال‌های ۱۹۹۰ ۲۰۲۰ و ۲۰۵۰ 
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17 شاخص بزرگ‌ترین لکه و LSI‏ شاخص JSS‏ سیمای سرزمین) در مقباس زیرحوضه‌ها در بازه زمانی 
۰ تا ۲۰۲۰ 


1 Number of Patch 

2 Patch Density 

3 Largest Patch Index 

4 Landscape Shape Index 


سال دوازدهم اثر تغییرات کاربری اراضی بر تولید آب در vee‏ ۳۹ 


نتایج ارزیابی تولید آب در حوضه آبخیز تراز با استفاده از اینوست در سطح پیکسل برای سال‌های 1۹۹۰ ۲۰۲۰ 
و ۳۰۵۰ در شکل ۵ یشان داده شده است. اسداللهى و همكاران' (۲۰۲۰) نيز ميزان توليد آب در حوضه آبخیز 
گر گانرود در استان گلستان را بين 55 تا ۵۳۳ ميلىمتر در هكتار برآورد كردند. مشخص كرديد که ميانكين توليد آب 
در حوضه در یک پیکسل ۰ متر برای سال‌های ۰ T°‏ ۰ به ترتيب AY VV‏ و AV‏ میلی‌متر (به ترتیب 
معادل ۸۵۰ 404 و 414 میلی‌متر در هکتار) است. كل حجم آب تولیدی قابل‌دسترس در منطقه در سال ۲۰۲۰ معادل 
است که به‌طور میانگین برابر با ۸۱۲ و ۱ مترمکعب در هکتار است. بدین ترتیب تولید آب در بازه زمانی ۰ نا 


هکتار) تغيبرات افرايشى داشته باشد, 


(mm) 


(mm) 


n 
” p” 
0 S 5 0 


شکل ۵- نقشه‌های تولید آب (mm/pixel)‏ منطقه مطالعاتی در ۱۹۹۰ ۲۰۲۰ و ۲۰۵۰ 
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ساخته شده ها کشاورزی 


شکل 1 میانگین تولید آب برای هر کاربری در سال‌های ۱۹۹۰ ۲۰۲۰ و ۲۰۵۰ 


میانگین تولید آب قابل‌دسترس در کاربری‌های جنگل» مرتع. کشاورزی و ساخته‌شده‌ها در سال ۲۰۲۰ به ترتیب 
۷۱ ۵ 2 و ۳۱۲۳ مترمکعب در هکتار و در سال ۱۹۹۰ به ترتیب ۱۱۳۵ ۷۲۶ 0۵ و ۲۳۹۲ مترمکعب در 
هکتار محاسبه شد. ميانگین تولید آب برای کاربری‌های مذکور به ترتیب ۰۱۰۵۷ ۱۲۸۳ ۳۸۵ و ۳۱۲۶ مترمکعب در 
هکتار برای سال ۲۰۵۰ پیش‌بینی شده است (شکل 1). بر اين اساس می‌توان نتيجه كرفت که در حال حاضر بیشترین 
ميزان تولید آب در منطقه در کاربری ساخته‌شده‌ها و در وهله بعد به ترتیب در مراتع» جنگل‌ها و اراضی کشاورزی 
اتفاق می‌افتد. احمدی ميرقائد و همکاران (۱۳۹7) نيز میزان تولید آب در کاربری‌های جنگل» مرتع» کشاورزی و 
ساخته‌شده‌ها در حوضه آبخیز قره‌سو در استان گلستان را به ترتیب LOOT ALY‏ ۹۵۳ و ۱۱۱۵ مترمکعب در هکتار 
برآورد کردند. تبخیر و تعرق زياد و قدرت حفظ و نگهداری آب در جنگل‌ها ازجمله دلایلی است که باعث کاهش 
آب تولیدی قابل‌دسترس در جنگل‌های منطقه نسبت به مراتع و ساخته‌شده‌ها شده است. همچنین, ميزان نفوذپذیری 
SLE‏ توپوگرافی و مکان جغرافیایی کاربری‌ها نیز در اين اختلاف اثرگذارند. SY‏ و همکاران (۲۰۱۷) نیز تأکید 
کردند که جنگل‌ها باعث کاهش تولید آب و ساخته‌شده‌ها باعث افزایش OT‏ می‌شوند. مقایسه تولید آب قابل‌دسترس 
در سال‌های ۱۹۹۰ ۲۰۲۰ و ۲۰۵۰ حاکی از افزايش OT‏ در کاربری‌های مختلف دارد که ازجمله علل آن می‌توان‌به 
تغيير شرايط محیط‌زیستی اکوسیستم‌ها و تغییر نفوذپذیری خاک‌ها در اثر تبدیل کاربری‌ها اشاره کرد. افزایش 
رواناب‌های سطحی و کاهش ذخیره آب‌های زیرزمینی در منطقه در طولانی‌مدت ازجمله پیامد اين فرآیند خواهد بود. 
پلنگ و همكاران' (۲۰۲۱) تأييد کردند که تبدیل کاربری‌های جنگل و مرتع به ساخته‌ضده‌ها Gel‏ افزایش 


1 Yang et al. 
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شكل ۷- طبقه‌بندی زیرحوضه‌های منطقه مطالعاتى از نظر توليد آب (mm/ha)‏ در ۸۱۹۹۰ ۲۰۲۰ و ۲۰۵۰ 


زیرحوضههای منطقه از نظر تولید آب در چهار طبقه خیلی کم mm/ha)‏ ۵۰۰ <( کم mm/ha)‏ ۵۰۰-۱۰۰۰ 
متوسط mm/ha)‏ ۱۰۰۰-۱۵۰۰) و زياد mmfha)‏ ۱۵۰۰ >( طبقه‌بندی شدند (شکل (V‏ نتایج نشان داد كه در سال 
۰ از کل مساحت منطقه ۱۲۳۹۲ ۱۳۳ ۵۸۹۰ و ۵۱1 هکتاربه ترتیب به طبقات با تولید آب خیلی کې کې 
متوسط و زياد احتصاص می‌یابد. اين مقادیر در سال ۱۹۹۰ به ترتیب ۱۲۹۵7 BAVA‏ 1۳۵۳ و ۳۱۳۷ هکتار است و 
در سال ۲۰۵۰ به ترتیب ۰۱۲۰۵۷ ۸۰۱۸ 14۳۵ و 1410 هکتار خواهد بود. شکل ۸ نيز روند تغییرات کاربری اراضی 
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شكل ۸- درصد تغييرات كاربرى اراضی و تولید آب در منطقه مطالعاتی در مقیاس زیر حوضه‌ها در بازه زمانی 
۰ تا ۲۰۵۰ 
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این نتایج نشان مىدهد كه بخش‌های شمال شرقى حوضه نسبت به ديكر مناطق از نظر توليد آب داراى بازدهى 
بهترى هستند و كمترين بازدهى نيز در قسمت‌های جنوبى منطقه مشاهده شده است. بارندگی زياد وجود دره‌های و 
كودالهاى عمیق نفوذيذيرى خاک و تراكم پوشش گیاهی مناسب در نیمه شمالی حوضه ازجمله عواملى هستند که 
در حفظ و ذخیره آب در اين پهنه‌ها اثرگذارند و باعث افزايش تولید آب در اين بخش‌ها نسبت به So‏ مناطق 
حوضه شده است. در سطح زیرحوضه‌ها نیز زیرحوضه‌های واقع در شمال حوضه از نظر تولید آب دارای بازدهی 
بهتری هستند درحالی‌که زیرحوضه‌های واقع در جنوب و جنوب غرب منطقه کمترین بازدهی به خود اختصاص 
داده‌اند. 

نتایج GWR‏ نیز در شکل ٩‏ نشان داده شده است. نتایج حاکی از آن است که تولید آب با سنجه‌های تعداد لکه 
تراکم لکه و شاحص شکل سیمای سرزمین رابطه مکانی معنادار و معکوس دارد درحالىكه رابطه آن با شاخص 
بزرگ‌ترین لکه یک رابطه معنادار و مستقیم است. در اين رابطه ضریب تبیین RI‏ بين تولید آب با سنجه‌های تعداد 
لکه تراکم لكه. شاخص شکل سیمای سرزمین و شاخص بزرگ‌ترین لکه به ترتیب ۰/۸۶ ۰/۸۳ ۰/۸۵ و ۰/۸6 تعيين 
شد. بر این اساس می‌توان كفت که ویژگی‌های ساختاری سیمای سرزمین و الگوی کاربری حاکم بر منطقه می‌تواند 
تعبین‌کننده تولید آب باشد. ولدگبریل و یشیتلا" (۲۰۲۱) تغییرات خدمات اکوسیستمی بر اثر تغییرات زمانی و مکانی 
کاربری اراضی را مورد تأیید قرار دادند. فام و همکاران" (۲۰۱۸) نیز تأیید کردند که تغییرات کاربری اراضى و 
ویژگی‌های آن باعث تغییر خدمات اکوسیستمی می‌شود. در واقع با افزايش تعداد تراکم و شکل لکه‌ها در یک قطعه 
از سرزمین خردشدگی و از هم گسستگی کاربری‌ها افزايش بيدا می‌کند كه در نتيجه آن شرایط اکولوژیک اکوسیستم‌ها 
تحت تأثير قرار می كيرد و در نتيجه اين فرآیند تغییر و کاهش خحدمات اکوسیستمی ازجمله تولید آب اتفاق می‌افند 
درحالی که با افزايش مساحت لکه‌ها (شاخص بزرگ‌ترین لکه) یک روند معکوس ایجاد می‌شود. 


1 Weldegebriel and 2 
2 Pham et al. 
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شاخص بزرگ‌ترین لكه و ‘LSI‏ شاخص شكل سيماى سرزمين) بر اساس GWR‏ 


۶- جمع‌بندی 

اين مطالعه نشان داد که تغييرات کاربری اراضی طی ۳۰ سال گذشته در حوضه آبخیز تراز در استان خوزستان 
شدید بوده است؛ به‌طوری‌که بیش از 1۵۰0 هکتار از اراضی جنگلی و مراتع منطقه تخریب و به اراضی کشاورزی و 
اهامای شدای Seats‏ مر كور سال اتف ود ۲۶۵۲ کک از رای کی .در 
صورت ادامه شرایط محیطزیستی کنونی حاکم بر حوضه کاسته شود. نتايج مطالعه تأييد کرد که میزان تولید آب در 
منطقه در کاربری‌های مختلف متفاوت است؛ به‌طوری که بیشترین و کمترین میزان تولید آب به ترتیب در کاربری‌های 
ساخته‌شده‌ها و کشاورزی برآورد شده است. همچنین ثابت شد که ویژگی‌های سیمای سرزمین در تغییر تولید آب 
اثرگذار است. در اين رابطه تولید آب با افزايش سنجه‌های تعداد. تراکم و شکل لکه‌های کاهش و با افزایش مساحت 
کاربری‌ها (شاخص بزرگ‌ترین لکه) افزایش می‌یابد. بر اين اساس می‌توان كفت که با افزايش توسعه کاربری‌ها 


ویژگی‌های سیماهای سرزمین تغيير و باعث افزايش گسستگی و خردشدگی زیستگاه‌ها می‌شود که در پی آن تغییر 
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کاربری‌ها و ارزيابى اثرات bol‏ در منطقه مطالعاتى بسيار حائز اهميت است. 


تشكر و قدردانى 
اين مقاله مستخرج از طرح يزوهشى يسادكترى به شماره 606 تحت حمايت صندوق حمايت از 
پژوهشگران و فناوران كشور است. نويسندكان از حمايت مسئولين و دستاندركاران صندوق و همجنين ويراستاران 
وداوران محترم به خاطر ارائه نظرات و دیدگاه‌های ارزشمند در راستاى افزايش كيفيت اين مقاله قدردانى شايسته به 
عمل می‌آورند. 
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